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Abstract 

We study the pure gauge model on a lattice manifold with triv- 
ial fundamental homotopy group, homotopically equivalent to an S4. 
Monopole loops may fluctuate freely on that lattice without restric- 
tions due to the boundary conditions. For the original Wilson action 
on the hypertorus there is an established two-state signal in energy 
distribution functions which disappears for the new geometry. Our 
finite size scaling analysis suggests stringent upper bounds on pos- 
sible discontinuities in the plaquette action. However, no consistent 
asymptotic finite size scaling behaviour is observed. 



1 Motivation 

Pure gauge lattice QED is a prototype model of gauge interactions of abelian 
nature and was among the first gauge systems studied in Monte Carlo sim- 
ulations |], 0, [3|]. It is also a starting points of the discussion of the more 
realistic lattice field theory containing dynamical fermion degrees of freedom. 
Unfortunately the status of nonperturbative studies of full QED is still con- 
troversial. In this situation it is important to clarify the fixed point properties 
i.e., the critical properties and critical exponents of the pure gauge system. 

From the early beginnings this field theoretic lattice model has posed 
many interesting questions about its continuum limit. The compact formu- 
lation of the pure gauge model has a confining strong coupling phase and a 
spin-wave (Coulombic) phase. The order of the phase transition has been 
controversial throughout the years. However, it was generally accepted that 
the formulation with the Wilson action on a finite lattice leads to a clear two- 
state signal at the transition, indicating a first order deconfinement phase 
transition Q. 

Adding extra terms to the action introduces new couplings. In this larger 
space the nature of the phase transition may change. Allowing for an adjoint 
coupling 7 (multiplying the plaquette variable in the adjoint representation 
[§]) a clear increase of the energy gap towards positive 7 was observed [|]. 
Based on a detailed study of discontinuities of the plaquette expectation value 
in the j3 — 7 plane and, using a power law scaling ansatz for the vanishing 
of discontinuities along the first order phase transition line, the location of a 
possible change of the order of the deconfinement phase transition to a second 
order phase transition line was estimated @, ||. From that extrapolation 
the infinite volume location of the so-called tricritical point in pure gauge 
compact QED was expected to lie at slightly negative values of the coupling 
7. Thus the theory on the Wilson line was thought to have a weak first 
order phase transition, however close to a tricritical point in the extended 
parameter space. 

Let us stress here, that up to now no direct simulation of pure compact 
QED did unanimously show the existence of a second order phase transi- 
tion point. It is well known from statistical mechanics systems that sizeable 
crossover effects from a discontinuous to a continuous behaviour are present 
on finite lattices in the vicinity of a tricritical point. Such effects are also 
present in pure gauge compact QED even at reasonably large negative val- 



ues of 7 and they spoil attempts to extract critical exponents of the system 
using finite size scaling methods. However, due to the central role of the 
model as a starting point of full QED and due to its prototype character in 
lattice gauge studies, it would be very important to know e. g. the value 
of the correlation length divergence exponent v. A possible determination 
then could tell us whether D = 4 pure gauge compact QED might be in the 
universality class of the free theory with v = 0.5, or not. We may remark 
here that MCRG studies of the model showed in the vicinity of the Wilson 
line a flow of couplings, which in spite of the two state signal resembled the 
flow diagram of a second order critical point quite well |7| . 

The nature of the deconfinement phase transition appears to be inherently 
connected to topological excitations present in this formulation (where the 
gauge fields assume their values in the group U(l), not in the algebra). For- 
mulated on the dual lattice these are closed loops of monopoles. In several 
investigations a strong correlation between the behaviour of the monopole 
loop density with the internal energy variable was observed ||, §, [|, [l(], [TT|1 . 
If the monopole loops are forbidden or suppressed by further terms in the 
action || [TTJ] a shifting of the phase transition towards smaller values of the 
plaquette coupling (3 was found. 

Here we concentrate on a modification of the simulation, which appears to 
be softer than introducing additional terms to the action, and which should 
become irrelevant in the thermodynamic limit. Instead of choosing periodic 
(hypertorus) boundary conditions we work on a lattice topology homotopic 
to the surface S 4 of a 5D hypersphere. In the infinite volume limit the two 
versions should lead to equivalent results. 

There is a subtle difference, which, however, may be important. The 
monopole loops are nonlocal objects. The Monte Carlo update algorithms are 
intrinsically local (except for scarce attempts towards multigrid updating). 
This is related to the gauge symmetry nature of the field variables. On the 
hypertorus closed loops lie in equivalence classes defined by the fundamental 
homotopy group of the torus: n x (T 4 ) = Z 4 . This puts some restrictions on 
the fluctuation behaviour of such loops. Hypercooling a hot configuration 
shows that monopole loops wrapped around the torus may got stuck like 
rubberbands [§. Recently this behaviour was studied and classified in more 
detail [E|. On the other hand, the fundamental homotopy group of the 



surface of a sphere is (for D > 2) trivial: n 1 (S' 4 ) = /. There monopole loops 
may contract and disappear without restriction. 



Preliminary results with half the statistics and for 7 = have been pub- 
lished earlier F 



2 Action and Lattice Geometry 

The action on the lattice A is given 

s = (3 £ Mu P ) + 7 E Mu 2 P ) (i) 

PgA PeA 

where Up denotes the usual plaquette variable and we consider all plaquettes 
with equal weight (see, however the discussion at the end). 

For the construction we choose the 4-manifold as the boundary of a 5D 
hypercube H 5 (L) with linear extend L (L 5 sites). We denote this boundary 
by SH^(L); it consists of all sites which have at least one cartesian coordinate 
value 1 or L. The manifold SH±(L) has 



sites 

links 

plaquettes 

3-cubes 

4-cubes 



10(L-l) 4 + 20(L-l) 2 + 2 

40(L-1) 4 + 40(L-1) 2 

60(L-1) 4 + 20(L-1) 2 (2) 

40(L - l) 4 

10(L-1) 4 



thus ratios of e.g. ^ii n ks/ n sites f° r the SH^L) lattice approach for large L 
asymptotically the values for the hypercubic geometry. Whereas the hyper- 
torus lattice T 4 (L) is self-dual, SH^L) is obviously not. For an identification 
of the monopole loops we did construct the dual manifold SH' A (L) (which has 
the same homotopy properties) as usual, identifying the dual sites with the 
centers of the hypercubes of the original lattice. This was done to check the 
dynamics of the updating algorithm. 

We define the volume of the lattice V = |A| as the number of plaquettes 
and we use 7< as the linear size variable; as finite size scaling is asymptotic 
this definition is as arbitrary as any other. In the computer programs the geo- 
metrical properties have been implemented by tables providing vectorization 
possibilities. 



3 Simulation and Coupling Parameters 

We performed long runs with an integrated statistics of 500K up to 1500K 
sweeps for each lattice size, using a combination of 3-hit Metropolis updating 
with one overrelaxation step. 

Lattices SH(L) were studied for L = 4, 6, 8, 10 having roughly the same 
number of variables as 6 4 . . . 16 4 lattices with torus geometry. Most of our 
data was taken for values of 7 = and /3 close to the transition. 

Simulations were done for various coupling values close to the pseudo crit- 
ical values. For given lattice size the histograms from different couplings were 
combined to one multi-histogram according to the technique of Ferrenberg 



and Swendsen [14|. Here we introduce statistical mechanics notations and 



define the partition function by 

Z L (l3) = Y,PL(E)e W (-[3E), (3) 

E 

where we denote the internal energy function by E = —J2pRe{Up). Note 
that E is extensive and negative for positive /5-values. Correspondingly we 
can introduce energy densities e = E/V, which up to a negative sign corre- 
spond to the plaquette operator. Here p(E) denotes the density of states for 
each value of E and is determined from the combined individual distribu- 
tions. In our simulation we determine expectation values of the specific heat 



and two cumulants KIR IT6 



CviP,L) = U{E-(E)f) 



V 

l ((E 2 -(£ ) 2 ) 2 ) 

3 (E 2 y 

{{E-(E)Y) 
((E-(E)W 



VboML) = -± "- ( ^7" , (4) 



Our definition of Vbcl differs by an additive constant | from that of [|i~5| 
The specific heat Cy develops a maximum, while Vbcl and U4, both develop 
minima as functions of /3 for given lattice size. We define pseudocritical 
couplings through the peak positions of these quantities and peak values of 
the considered quantities by their values at pseudocriticality. 

We observed quite long autocorrelation lengths at the corresponding pseu- 
docritical points for each lattice size, with peak values of Ti nt ranging from 



15 (for L = 4) up to 750 (for L = 10). This was one of the reasons for 
doubling the statistics and introducing overrelaxations steps as compared to 
the preliminary results presented in W3 . 



4 Results and scaling analysis 

4.1 Pseudocritical couplings 

Table 1 gives the pseudocritical couplings and the values of the cumulants at 
the peak positions. 

In fig. 1 we compare the multihistograms resulting from all data for 
given L, that is /i.l(/?, e = E/V) = Pl(E) exp (— j3E) at the values of the 
peak positions of Cy. These are the effective probability distributions for the 
action at these couplings. We find no indication of double peaks, in distinct 
contrast to results for the same action but standard hypertorus geometry 

We then performed runs at 7 = 0.2 where torus-results had identified even 
more prominent two-state signals. An analogous multihistogram analysis 
there showed pronounced two-state signals for the theory on SH lattices 
even on the smallest lattice sizes. Whenever this situation occurs, tunneling 
becomes scarce on larger lattices and the multihistogram analysis becomes 
impracticable. In this case one should rely on multicanonical simulations 
For that reason we did not attempt simulations for SH(L > 6) at this 



value of 7. 

Returning to 7 = we display in fig. 2 as a function of the variable 1/V a 
comparison of pseudocritical j3- values on SH lattices, namely /3(Cv), (3(Vbcl) 
and /3(U4), with results from simulations on the standard hypertorus geome- 



try for (3(Cv) [0- We find a strikingly different finite size scaling behaviour 
when comparing both sets of data. The pseudocritical /3-values of the SH 
lattices approach the infinite volume critical /3 c -value from above, while data 
on the hypertorus do so from below. Furthermore our data indicate in case of 
/3(Vbcl) a non-monotonic behaviour on the SH lattice, while on the standard 
hypertorus the approach seems to be monotonic for all quantities. 

In principle finite size scaling theory gives the asymptotic approach of the 
considered pseudocritical j3- values. If the deconfinement phase transition is a 
continuous phase transition the homogeneous free energy density fhom(P) — 



fhomVL-) scales with the scaling variable; with the usual definition of the 
correlation length divergence exponent v one has for large volume systems 
the asymptotic behaviour 

/3 cy i=f3 c + aV-^. (5) 

Alternatively, if the phase transition is of discontinuous nature (first order), 
we expect an expansion of the pseudocritical /3- values in powers of the inverse 



volume 1/V pO, |T 



f3 cy = (3 c + v + 0(l/V 2 ). (6) 

In fig. 2 we display a fit to the /3(CV)-values on the standard hypertorus 
according to the scaling law (g) including 1/V 2 corrections (dashed curve in 
the figure). This fit gives a consistent description of the data and results 
in the infinite volume /3 c -estimate (3 C = 1.01102(9). In the current context 
this demonstrates the consistency of the data with a finite size scaling be- 
haviour of a discontinuous transition in the considered volume size interval. 
As already mentioned, this was anticipated for the case of the hypertorus 
geometry. That fit alone is not sufficient, though. A definite statement 
on whether the deconfinement phase transition is of first order type would 
require the calculation of its latent heat and interface tension in the thermo- 
dynamic limit. Such an analysis is beyond the scope of our present work. 
Later on we argue, indeed, that evidence for a weaker first order or even 
continuous phase transition can be obtained from SH lattices. 

Returning to the theory on SH lattices we note that by construction the 
theory is not homogeneous. On the surface of the 5D hypercube there are 
inhomogeneities related to its edges i.e., the theory on the SH lattice exhibits 
less symmetries than formulated on a hypertorus. In general one then ex- 
pects additional finite size scaling violations which may be parametrized by 
a contribution fi n hom{P) to the free energy density on the SH lattice. In the 
thermodynamic limit fmhomift) will vanish. 

In absence of a rigorous theory we may estimate the dimensionality of the 
additional terms by considering quantities sensitive to the inhomogeneity of 
the lattice. For instance we might consider ratios of the number of plaquettes 
to the number of links (fj), or the number of links which are contained in 4 or 5 
different plaquette variables as compared to the standard case where each link 
is contained in 6 different plaquette variables. These considerations suggest 
a dimensionality of the inhomogeneous free energy density part oc l/Vz. 

7 



The usual ansatz for the partition function (the indices o and d denote or- 
dered and disordered phase) of a discontinuous system, neglecting interfacial 
effects, is 

Z(f3) = exp (-Vf d ) + Aexp (-Vf ). (7) 

Here the quantity A parametrizes the asymmetry. The free energy densities 
fo(fl) and fd(P) can be expanded in powers of Af3 = (3 — (3 C 

Jo Jo,hom i Jo,inhom 

= c + eoAp + V-^rio + aoA^+OttAP) 2 ) (8) 

(and correspondingly for /<$). In the infinite volume limit f = fd = c at 
A/3 = 0. At finite volume we find at the metastability point, where both 
contributions have same weight, 

\nA + e A(3 + V~^ (r] + a A(3) = 

e d Ap + V-^Vd + ffdA/3) + 0((A/?) 2 ) . (9) 

Solving for A(3 we then expect the finite size scaling law of (H) to be replaced 
by the form 

Pc,V = Pc+^T + ^ + 0(V^), (10) 

V 2 V 

if the phase transition was of discontinuous nature. It turns out that the 
dominant term now is 0(1 /Vz). 

In fig. 2 we show fits of the pseudocritical /3-values for the SH lattices 
according to ( |T0f) (dotted curves). These fits result in a positive contribu- 
tion proportional to 1/Va and a negative contribution proportional to 1/V, 
which actually is of the same order of magnitude as in case of the hyper- 
torus. Thus the observed non-monotonic behaviour observed in the quantity 
(3(Vbcl) could be explained by the interplay of a contribution related to the 
inhomogeneity, entering with the opposite sign, as compared to a contribu- 
tion ex 1/V, which may be attributed to the homogeneous part of the free 
energy density. The obtained infinite volume /3 c -estimates from these fits 
(3 C (CV) = 1.01133(13) , (3 c (V B cl) = 1.01134(15) and f3 c (UA) = 1.01125(16) 
are consistent with the result from the hypertorus, though a tendency to- 
wards values slightly higher by about one standard deviation is noticeable. 
Again we emphasize that this analysis is merely a statement on the consis- 
tency of different finite size scaling fits in the considered and limited volume 



size interval i.e., both regularisations on the hypertorus and the SH lattice 
support a unique infinite volume /5 c -value. As far as the SH lattice is con- 
cerned, we conclude, that the behaviour of the pseudocritical /3-values is 
strongly affected by contributions stemming from the inhomogeneity of the 
system. Consequently we have to pay special attention to the effect of these 
contributions when in the later analysis we study finite volume estimators of 
energy gaps at the deconfinement phase transition. 

We have also analyzed the SH lattice pseudocritical /3-values assuming 
the finite size scaling form (g) corresponding to a critical point. Including 
data with L > 8 into the fit we obtain ^-values of about v ~ 0.53(20) and 
/^-values consistent with the above determinations. However, as discussed 
such a value of v close to 0.5 can just be mimicked by contributions related to 
the inhomogeneity of the system fllPl). We therefore cannot trust these fits, 
as long as these contributions are sizable and likely to be present. This obser- 
vation also indicates that an unambiguous extraction of critical exponents of 
the deconfinement phase transition, if it were in fact a continuous transition, 
would require much larger lattice sizes, or alternatively, a formulation of the 
theory with smaller finite size scaling violations. 

4.2 Cumulant peak values 

The peak values of the specific heat and cumulants defined in @ give in 
the thermodynamic limit unambiguous signals on the order of the considered 
phase transition. If the phase transition is of discontinuous nature one finds 
in the thermodynamic limit L — > oo on general grounds 



V 4 

( e 2 _ e 2)2 
V BC L,min = ~ 12 ° (eo -p , (11) 

U 4,min ~ ~ *■ • 

Here e and e^ denote the discontinuous values of the energy density at the 
critical point. In the context of the present paper — e corresponds to the 
plaquette expectation value in the Coulomb phase, — e^ to its expectation 
value in the confinement phase. For a continuous phase transition e = e^. 
In this case U± m i n can be different from 1. 



Again finite size scaling theory predicts for discontinuous and homoge- 
neous systems, that finite size corrections to above quantities on finite vol- 
umes follow asymptotically an expansion in powers of the inverse volume 
1/V, similar to the expansion written down in (^). In fig. 3a we display a 
comparison of results on the hypertorus (triangles) with results on SH lat- 
tices (circles) for the quantity Cy jmax (L) JV as a function of the parameter 
1/V. This quantity can be directly interpreted as a finite volume estimator 
of a possible infinite volume discontinuity in the internal energy or the pla- 
quette expectation value, correspondingly. We observe for both geometries 
a deviation from the leading asymptotic finite size scaling with terms pro- 
portional 1/V i.e., higher order terms in the corresponding expansions are 
important. It is therefore rather doubtful, whether a fit to the data, even 
including higher orders, could lead to a unambiguous determination of the 
infinite volume discontinuities. 

Here we are not that ambitious. We note, however, that the SH data 
exhibit much lower values for the finite volume estimator of the discontinuity 
on comparable volume sizes. Thus, if the observed finite size scaling on the 
SH lattice stays monotonic in the variable 1/V up to to the thermodynamic 
limit, then the SH results predict rather stringent upper bounds on the values 
of possible discontinuities at the phase transition point. In light of the above 
discussion on possible effects related to the inhomogeneous part of the free 
energy density we have to be rather careful. We therefore performed in 
our analysis several fits to the observable Cv, ma x(L)/V including 1/V', 1/V 
and also corrections proportional to 1/V 2. In none of these cases we found 
indication of a non- monotonic behaviour. The dotted curve of fig. 3a displays 
as an example such a fit with the form 

Cv.max __ (e — Gd) a b 

Thus we convert the value obtained on the SH( 10)-lattice, Cy tTnax {lQ)/V = 
.0000528(13), to an upper bound on the possible discontinuity in e d — e in 
the thermodynamic limit 

e d - e < 0.0145(2), (13) 

which then should be satisfied by pure gauge compact QED with the Wilson 
action in the hypertorus geometry as well as on SH lattices. We remark here 
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that such a bound is more stringent than any gap observed or extrapolated 
from torus results 0; simulations on the hypertorus on 16 4 lattices typically 
result in estimates of the discontinuity of about 0.030 on that lattice size. 
In a study of the discontinuity in the j3 — 7 plane of couplings a power law 
ansatz suggested an infinite volume value of the above discontinuity of about 
0.016 0. However the statement here is that (|T3| ) constitutes a bound and 
that the true value of the discontinuity is likely to be still smaller. Thus 
compact pure gauge QED on SH lattices provides evidence that the first 
order deconfinement phase transition on the Wilson line is at least weaker 
than previously thought. 

In fig. 3b we also display a comparison of results on the hypertorus 
(triangles) with results on SH lattices (circles) for the quantity Vbcl,™,™- 
The remarks concerning deviations from asymptotic finite size scaling in the 
preceding paragraph are appropriate here as well. Performing various fits 
of the already mentioned forms (dotted curve in fig. 3b) we do not find a 
non-monotonic behaviour. Again the observed monotonic behaviour of the 
quantity VBCL,min lead to an upper bound on possible discontinuities. On 
the SH (10) lattice (cf. table 1) we obtain the inequality 

e d -e < -0.022(1) e d . (14) 

The plaquette expectation value in the disordered state —e d can be easily 
bounded from above; numerical simulation results give — e d < 0.66. Inserting 
this number we again obtain e d — e < 0.0145(2), which is precisely the bound 
derived in (|13|). Finally we display in fig. 3c the quantity U4. It does not 



exhibit asymptotic finite size scaling behaviour for either version of the theory 
and it is not possible to decide whether it approaches a value of 1 (for a first 
order phase transition) or larger (continuous transition). 

All together our findings indicate that asymptotic finite size scaling for 
both formulations, on the hypertorus and on the SH lattice, is very far away 
from the lattice sizes, which we are able to investigate in our present studies. 
This is not surprising for the theory on the hypertorus as we know from 
studies in much simpler models (like Potts models) with first order phase 
transitions, that asymptotic behaviour may set in only on very large lattices 



2^j and in general little is known about the onset of finite size scaling for 
gauge theories. On SH lattices additional finite size corrections are caused 
by the inhomogeneity. 
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5 Discussion 

The studied geometry - the boundary of a 5D hypercube - has less symmetry 
than the hypertorus. Asymptotically it amounts essentially to a sum over a 
set of ten hypercubes (like in 2D it corresponds to the six faces of the cube, 
the boundary contributions become negligible). On finite lattices the reduced 
translational and lattice rotational symmetry gives rise to large additional 
contributions of finite size scaling violating operators. Thus the asymptotic 
scaling behaviour sets in only on very large lattice sizes. A possible cure 
for the observed finite size scaling violations may be the introduction of 
weight factors i.e., considering the metric for a truly curved topology like the 
S4. Studies in this direction including also terms proportional to 7 into the 
action are under progress and we expect from those studies clarification on 
the critical properties of pure gauge compact QED close to the Wilson line. 
It remains a remarkable fact that the clear two-state signal of earlier 
Monte Carlo simulations disappears when one changes the lattice geometry 
with respect to their properties at the boundary. This finding corresponds 
to a rather small bound on possible discontinuities in the plaquette operator, 
which, based on the observed monotonic behaviour of finite volume discon- 
tinuities, is suggested by the present analysis. It also remains a remarkable 
fact, that the introduced soft modification of the theory leads in compari- 
son to the results for the theory on a hypertorus to such large effects. This 
may be related to the influence of closed loops for the dynamics of the phase 
transition and may be a special property of gauge theories containing the 
compact U(l) gauge group. One may wonder whether this observation can 
be of relevance for other systems with loop excitations, including systems 
with fermionic fields. 

Acknowledgement : 

In the final stage of our work we profited from discussions with Jiri Jersak. 
C.B.L. wants to thank Martin Liischer for a discussion. T.N. is indebted to 
Thomas Lippert and Klaus Schilling for a useful discussion on the subject. 



12 



References 

[1] M. Creutz, L. Jacobs and C. Rebbi, Phys. Rev. D 20 (1979) 1915. 

[2] B. Lautrup and M. Nauenberg, Phys. Lett. 95B (1980) 63. 

[3] T. A. DeGrand and D. Toussaint, Phys. Rev. D22 (1980) 2478. 

[4] J. Jersak, T. Neuhaus and P. M. Zerwas, Phys. Lett. 133B (1983) 103. 

[5] G. Bhanot and M. Creutz, Phys. Rev. D 24 (1981) 3212. 

[6] H. G. Evertz, J. Jersak, T. Neuhaus and P. M. Zerwas, Nucl. Phys. B 
251 [FS13] (1985) 279. 

[7] C. B. Lang, Phys. Rev. Lett. 57 (1986) 1828; Nucl. Phys. B 280 (1987) 
255. 

[8] V. Grdsch et al, Phys. Lett. 162B (1985) 171. 

[9] J. S. Barber, R. E. Shrock and R. Schrader, Phys. Lett. 152B (1985) 
221; J. S. Barber and R. E. Shrock, Nucl. Phys. B 257 [FS14] (1985) 
515. 

[10] G. Bhanot, T. Lippert, K. Schilling and P. Ueberholz, Int. J. Mod. Phys. 
C 3 (1992) 1109; Nucl. Phys. B (Proc. Suppl) 30 (1993) 912. A. Bode, 
T. Lippert and K. Schilling, Nucl. Phys. B (Proc. Suppl.) 34 (1994) 
549. 

[11] V. G. Bornyakov, V. K. Mitrjushkin and M. Miiller-Preussker, Nucl. 
Phys. B (Proc. Suppl.) 30 (1993) 587; V. G. Bornyakov, V. K. Mitr- 
jushkin, M. Miiller-Preussker, and F. Pahl, Dirac sheets and gauge fixing 
in U(l) lattice gauge theory, preprint HU Berlin-IEP-93/2, 1993. 

[12] W. Kerler, C. Rebbi and A. Weber, Phase structure and monopoles in 
£7(1) gauge theory, BU-HEP 94-7, 1994. 

[13] C. B. Lang and T. Neuhaus, Nucl. Phys. B (Proc. Suppl.) 34 (1994) 
543. 



13 



[14] A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 61 (1988) 2635; 
Phys. Rev. Lett. 63 (1989) 1195; Computers in Physics Sept. /Oct. 
(1989) 101. 

[15] M. S. S. Challa, D. P. Landau and K. Binder, Phys. Rev. B 34 (1986) 
1841. 

[16] K. Binder, Finite size effects at phase transitions, in Computational 
Methods in Field Theory, Lecture Notes in Physics 409, edited by H. 
Gausterer and C. B. Lang, page 59, Springer- Verlag, Berlin, Heidelberg, 
1992. 

[17] V. Azcoiti, G. D. Carlo and A. F. Grillo, Phys. Lett. 238B (1990) 355. 

[18] B. Berg and T. Neuhaus, Phys. Lett. B 267 (1991) 249; Phys. Rev. 
Lett. 68 (1992) 9. 

[19] T. Neuhaus, Multicanonical simulation of compact QED, in preparation. 

[20] C. Borgs and R. Kotecky, J. Stat. Phys. 61 (1990) 79. C. Borgs, R. 
Kotecky, and S. Miracle-Sole, J. Stat. Phys. 62 (1991) 529. 

[21] J. Lee and J. M. Kosterlitz, Phys. Rev. B43 (1991) 3265. 

[22] A. Billoire, T. Neuhaus, and B. Berg, Nucl. Phys. B 396 (1993) 779. 



14 



Tables 



Table 1: Peak positions and values of the measured cumulants. 



L 


^V,max 


VBCL,min X 10" 


U 4,min 


4 


2.21(1) 


-1.631(7) 


2.739(8) 


6 
8 
10 


4.96(13) 

10.08(32) 

20.91(62) 


-0.45(11) 

-0.23(8) 

-0.17(5) 


2.625(30) 
2.455(43) 
2.263(37) 




L 


P(Cv) 


P(Vbcl) 


P(U*) 


4 
6 
8 
10 


1.01536(22) 
1.01370(16) 
1.01263(6) 
1.01216(1) 


1.01194(29) 
1.01326(18) 
1.01257(6) 
1.01215(1) 


1.01688(37) 
1.01401(19) 
1.01273(7) 
1.01218(2) 
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Figure 1: The multihistograms /i.l(/3(CV), e) resulting from all data for given 
L at the values of the pseudocritical points as defined by the peak positions 
of the specific heat. Shown are data on SH (4) , SH (6) , SH (8) and SJf(lO) 
lattices. 
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Figure 2: The pseudocritical /3-values as a function of the parameter 1/V. 
We compare results for (3{CV) (circles), /3{Vbcl) (triangles) and /3(L r 4) 
(crosses) on the SH lattice with results from simulations for /3(CV) (triangles 
down) on the hypertorus. The data on the hypertorus have been obtained 
in multicanonical simulations [19j on lattices ranging in volume between 6 4 
and 12 4 . On the hypertorus we include also a data point from |10| on the 16 4 
lattice. The dashed curve represents a fit according to (|5|), the dotted curves 



represent fits according to (|10|). The full symbols represent infinite volume 
extrapolations of f3 c . 
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Figure 3: The peak values of (a) specific heat Cy and of (b) the cumulants 
Vbcl and (c) [/ 4 vs. 1/V. We compare data on the hypertorus (triangles) 
with data on SH lattices (circles). The dotted curves are explained in the 

text. 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407005vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407005vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9407005vl 



In h, (fi(CV) , e) 



' i^T*^ 






i ' i ' i 




L8 


\ 6 


\ 4 


- 


If 10 \ 






\ 


/ J 


i , i 


II 


1 i ,1 


i , i , i \ 



-.68 -.66 -.64 



-.62 



-.6 



-.58 -.56 

e 



-.54 



FIG. 1 



Figure 1: The multihistograms /ii(/3((7y), e) resulting from all data for given 
L at the values of the pseudocritical points as defined by the peak positions 
of the specific heat. Shown are data on SH (A) , SH (6) , SH (8) and Sfl^fO) 
lattices. 
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Figure 2: The pseudo critical f3- values as a function of the parameter 1/V. We 
compare results for f3(CV) (circles), /3(Vbcl) (triangles) and (3(U4:) (crosses) 
on the SH lattice with results from simulations for fi(CV) (triangles down) 
on the hypertorus. The data on the hypertorus have been obtained in mul- 
ticanonical simulations [19] on lattices ranging in volume between 6 4 and 
12 4 . On the hypertorus we include also a data point from [10] on the 16 4 
lattice. The dashed curve represents a fit according to (6), the dotted curves 
represent fits according to (10). The full symbols represent infinite volume 
extrapolations of f3 c . 
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Figure 3: The peak values of (a) specific heat Cy and of (b) the cumulants 
Vbcl and (c) U^ vs. 1/V. We compare data on the hypertorus (triangles) 
with data on SH lattices (circles). The dotted curves are explained in the 
text. 
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